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Abstract 

We study a tensorial exponential transformation of a three-dimensional metric of space-like 
hypersurfaces embedded in a four-dimensional space-time, 7^- = e! Lihm6m e^ k e~ eihn6n , where 
4>k are logarithms of the eigenvalues of 7^, 9k are rotation angles, and ey/. is a fully anti- 
symmetric symbol. Evolution part of Einstein's equations, formulated in terms of 4>k and 
9k, describes time evolution of the metric at every point of a hyper-surface as a continuous 
stretch and rotation of a local coordinate system in a tangential space. The exponential 
stretch-rotation (ESR) transformation generalizes particular exponential transformations 
used previously in cases of spatial symmetry. The ESR 3+1 formulation of Einstein's 
equations may have certain advantages for long-term stable integration of these equations. 
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1 Introduction 



3+1 formulations of general relativity (GR), commonly used in numerical integration, break 
Einstein's equations onto constraints that must be solved to obtain initial conditions on an 
initial space-like three-dimensional hypersurface, and a Cauchy part, an evolution problem 
that is then integrated forward in time 121- It is well known that integration in time is 
often unstable and terminates prematurely. In particular, this happens frequently in numerical 
calculations of black hole collisions which are expected to be a major source of gravitational 
radiation gravitational observatories [3] • There is no complete understanding of the causes of 
instabilities. Gauge and constraints instabilities of Einstein's equations are definitely a part of 
the problem 4 . Bad choice of a numerical scheme may be another. Einstein's equations are 
highly non-linear, and finding a numerical method capable of stable long-term integration of 
these equations remains an outstanding problem in numerical GR. 

In , we gave an example of an asymptotic non- linear instability which arises in a simple 
scalar hyperbolic equation, g tt = g xx — | {agf + (3g x + ^gxdt), whose non-linear term mimics a 
quadratic non-linearity of much more complex Einstein's equations. This equation is well-posed 
and it does not suffer from gauge or constraint instabilities. Stable numerical schemes applied 
to this equation converge to exact solutions at any fixed time t when numerical resolution is 
increased. However, numerical integration can become asymptotically unstable if the resolution 
is kept fixed and t is increased, t — > oo. The instability is caused by breaking down of a delicate 
balance between linear and non-linear terms in the discretized versions of the equation. An 
important point is that the breakdown takes place without violating second-order accuracy of 
a discretization. Therefore, having both convergence at finite t and an asymptotic instability 
at t — > oo is not a contradiction. 

We studied in an exponential transformation g = e^, which maps < g < oo onto 
— oo < <j> < oo and leads to a modified equation <p tt = <fi xx — {a + ^)4> x — {13 — 1)4> X — j(p x 4>t, 
for a logarithmic variable <f>. We demonstrated that using the modified equation instead of the 
original one dramatically improves the accuracy and long-term stability of numerical integra- 
tion. 

There is no reason to believe that asymptotic instability is limited to scalar equations 
and that it cannot occur in more complex systems such as equations of GR. We want to 
investigate whether discretized Einstein's equations are prone to such an instability, and if 
numerical integration of these equations can be improved and prolonged by an exponential 
transformation of variables. As a first step in this direction, we show that it is possible to 
formulate a tensorial exponential transformation for a three-dimensional metric 7^ of space- 
like hypersurfaces embedded in a four-dimensional space-time (Section 2). The transformation 
leads to a 3+1 formulation of Einstein's equations which describes the evolution of jij at 
every point of a hypersurface as streth and rotation of a local cartesian coordinate system in a 
tangential space (Section 3). 
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2 Exponential transformation of a three-dimensional metric 



2.1 Metric in terms of stretch and rotation variables 

To carry out a tensorial exponential transformation, we write a three-dimensional metric of a 
space-like hyper surface, dh 2 = ^ijdx l dx^ , as 

H3= A \k D kiAi 3 , (2.1) 

where D{j = 5ijX{ is a diagonal matrix of eigenvalues of jij, and Aij is the orthogonal matrix of 
rotations, A]- = Aj, 1 ; superscript f denotes a matrix transposition, A]- = An, and 4 m ji m , = 
5ij. Decomposition (|2. 1|) is always possible for a symmetric matrix. 

Since the three-dimensional metric is positive definite, < Aj < oo, we can introduce 
logarithms of eigenvalues 4>i = In Aj and write 

Dij = Sije^, -oo < 0i, 02, 03 < oo. (2.2) 

The transformation maps eigenvalues Aj onto a (— oo,oo) interval; dh 2 = corresponds to 
4>k ->■ -oo. 

Rotation matrix can be further expressed as exponentiation of an anti-symmetric matrix. 
We write 

A^ = exp(My), Mij = e ijk 9 k , -oo > 9%, 9 2 , 9 3 < oo, (2.3) 

where 9 k is a rotation vector dual to Mij, and eyfc is a fully anti-symmetric symbol (matrices 
= e *j fc > ^ = 1 5 2, 3, are generators of a rotation group). Using new variables, we can rewrite 
(J23J) as 

7y = e tikm6m e $k e -ejkn6n_ ( 2 .4) 

Equation (|2. 1|) or (|2.4|) describes 7y as a result of stretch and rotation of a cartesian coordinate 
system with a metric 5ij carried out in a tangential space at every point of a hypersurface. 
Transformation (|2.1[) is unique if all three eigenvalues of jij are different. If two or three 
eigenvalues are degenerate, then one of the angles of rotation or all three of them cannot be 
uniquely defined. However, quantities which have physical meaning, such as angular velocities 
and acceleration, remain meaningful. The evolution equations for 0j and 9{ which we formulate 
below are unique (we address this in detail at the end of the paper). In what follows, we will 
refer to <pi, B%, and their certain derivatives as exponential stretch rotation (ESR) variables. 
Our aim will be to rewrite Einstein's equations in terms of ESR variables. 

Before proceeding any further, we must discuss index notation and summation rules used 
in the paper. These rules are different from these commonly used in GR because ESR variables 
and matrices A^ do not transform like vectors or tensors under curvilinear coordinate trans- 
formations, and it usually takes a combination of several ESR variables to form a covariant or 
a contravariant object. In index notations, we will not distinguish between upper and lower 
indices in ESR variables and their derivatives. For example, A l j will be a matrix with the same 
elements as A^ or A lJ . Summation in formulas involving ESR variables is implicitly assumed 
over a repeating index. An index in a formula may repeat more than two times. However, 
there will be one (and only one) very important exception from this rule: summation must 
never be carried over a repeating index if the index is present only once on one side of an 
equation and more than once on the other side. Equation Q2.2JI above is an example of the 
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summation exception rule. The is no summation over i in this formula. Another example will 
be an equation a* = b] m c k idVj kl , where summation must be carried out over indices m, k, and 
I, but no summation is carried over i. However, in the equation a\ = b l lm c kl d™ kl summation 
must be carried over m, k, /, and i on the right-hand side, and over i on the left-hand side. 
These modified rules let one to simplify and unambiguously perform manipulations with equa- 
tions involving ESR variables, and to express tensorial objects in terms of ESR variables in a 
compact and unambiguous way. 

In what follows, it will be convenient to use auxiliary variables 



m = j, o = yJd\ + el + el -k%<i. (2.5) 

Since e ni ke n j m = 5ij5km - $imfaj, we have 

Mfj = e ink e k e njm 6 m = 9 2 (rnVj ~ <%) ■ (2-6) 
Rodrigues formula then gives the rotation matrix Aij in terms of 9 k , 

, , sin 9 , , 9 1 — cos 9 



5ij cos 9 + eijkrjk sin 9 + rjirjj (1 — cos 9) . 



Using (|2.7[) . we can write a closed expression for the metric (|2.1[) as 



7ij =5 ij e <pt cos 2 9 + ( e 9i - j e ijk rj k sin 9 cos 6*+ 

einkejnmVmVke*" sin 2 + (V 1 + e*A r/^ cos (9(1 - cos 9)- ( 2 -^ 
(^n^j + ejnkVi) e 4 ' n r] k r] n sm9(l - cos9) + e^rj^mq^l - cos9) 2 . 



Since 



'fij — A m {€^' A m j , 7'°' — (-^mi^ A m j^ — A rn j i A m je ^ , (2-9) 
formula Q2.8J1 with (pi — * —(pi gives a closed expression for ^ . 

2.2 Differentials of the metric 

We now must derive the formulas which relate first and second differentials of the metric, d'jij 
and d 2, jij, with first and second differentials of ESR variables, dcpi, d9i, d 2 (pi, d 2 9i. Differentia- 
tion of (|2.1j) gives 

d'jij — dA^ n D nm A m j -\- A, m D nm dA m j + A^ n dD nm A m j 
B k - n ie^ n A n j + A n ie^ n B k . n j ) d9 k + (A k ie^ k A k j 



where we introduced three new matrices B k which are first derivatives of a rotation matrix 
with respect to 9 k , 

dAij sin 9 , „ s 1 — cos I 



Sill (7 . . 1 - CUS (7 

^Ajij =1m _ = e *i fe 2T - + WkVj + OjkVi) 2 h 



r] k - s'm95ij + e ijn r} n [cos 9 — + rj^j sin 9 - 2 



sin \ / . „ „ 1 — cos v x x 
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In (|2,lll) and throughout the paper, symbol ";" does not mean a covariant differentiation, but 
simply used to separate indices of different nature. 

To invert (j2.1(Jj) and obtain d<p k and d6 k in terms of djij, we must rotate (j'i.lOj) into a 
coordinate system where 7^ is diagonal by multiplying l|2.10|) with Aij and Al- . The result is 

AindjnmA^ = C ki ijd6 k + <%e^d<^, (2.12) 

where 

Ck;ij = AinBj : . nrn D rn j -\- Di n B k - nm A m j — Ai n B k j n e^ >: ' -\- Aj^B^-ine^ 1 . (2.13) 
From orthogonality of A{j it follows that combinations Ai n B^ k . n - are anti-symmetric, 

-^m-^j.. n j — Ai n Bk-j n — Bfc-i n Aj n — Aj n B £, n j- (2-14) 

Therefore, C k -ij are symmetric trace-free matrices with all diagonal elements identically equal 
to zero, 

Ck;ij = Cjfcyi) Cfc ; ii = Cfc ; 22 = Cfc;33 = 0, (2-15) 

and because of this property dcj) k and d0 k in (|2.12() do not mix. We can use diagonal part of 
(l2~T2l to find 

dcf)i = e"^ A im A in dj mn . (2.16) 
Off-diagonal part of (|2.12j) gives a system of three linear equations 

Ai n d'y nrn Aj rn = Ck-^jdOk, i 7^ j, (2.17) 

which can be solved to find d9, L , 

d0i = C^j \€jmn\-A mr A n j : d'y r k, (2-18) 

where 

Cij — I Cinm I Cjr'jnm,) (2.19) 

and | . | is a module operation. Using matrices Ck-ij we can write 

djij = AkiAkje^dfik + A ni A mj C k]nm d6 k . (2.20) 
To derive expressions for d 2 (j)i and d 2 6i we differentiate ()2.12(l and obtain 



(2.21) 



where 



5 ij e^d 2 (t> i + C k .i j d 2 9 k = 

H-ij;nm = ~&ij&irn,5jn^ t i (2.22) 

-ryl _ ttI \ ttI n_/2 _ tt2 < tt2 icy nn\ 

rl -ij;nm ij',nm ' ji\nmi r ^ij;nm -"ij'jnm ' ■ n ji;nmi y^.^oj 

■^ij;nk — -l^n;im-^rmC k -, r j 6 '* (yB k -i m B n -j m + E kn -i m Aj m ) , (2.24) 
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Hfj-M = B k , im A jm (e^S jn - e^Sin) , (2.25) 



and we denoted second derivatives of a rotation matrix as 

_ d 2 Aij fss l-cos# sin 6 



cos a sin ( 
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tijrVk + eij n r] n (5kr - 3r] r r]k)) 



I sin \ 

eijnVrVkVn sin 9 - 5ijr] k r] r I cos 9 — j + (2.26) 

(sin $ 1 — cos \ 
~9 2 — ff 2 — / + 

/ sin 9 ( sin 6* 1 — cos ( 
mVjVkVr I cos 9 4 — 2- 
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Note, that H%. nm d<j) n d^> m make a non-zero contribution only to a diagonal part of (|2.21|) . 
whereas 'H 2 j- nm d(j)nd9m contributes only to an off-diagonal part. 
Diagonal part of (|2.21j) then gives us 

d 2 fr = e-+* (A lk d 2 lkn A m + H} l . nm d9 n d9 m ) - (d&) 2 . (2.27) 



Off-diagonal part of (|2.21j) gives 

d 9% = C^j \€j mn \ (A mr A nk d ■jrk + TC^j. nm d9 n d9 m -\-T~C^j. nm dcj) n d9 m ) , (2.28) 

where dj has been defined in (t2~TTIl) . The above formulas, (t2~TCl) . (I^TTHl , ^TF!\ . and 

relate first and second differentials of new variables <j) k and to first and second differentials 
of jij. There is a linear relation between d9 k , d(f) k and d'jij ( (|2.12|l ). Both ci 2 and (f 2 ^ depend 
linearly on second differentials d 2 jij and in addition involve squares of first differentials djij 
(see {Eg} ). 

2.3 Cristoffels and Ricci tensor 

We use 1)2.10(1 to write first derivatives of 7^ with respect to spatial coordinates as 

- A A + e^° (B A ■ + A B .\®Hm±. (2 29) 

and Cristoffel symbols 
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as 



r ij - 2 Aak i^ ai ~j^j + Aaj ~d^ ^ <PaAamAbiAb ^'^^ + 

-"mjai ~ c ■ n an- n bi 1J m\bn J Q^,j ' 

J -''m;aj ~ c ■ n an- n bj- LJ rn;bn J q i 

d8 m \ 



(2.31) 



From (|2.31|) we get contracted Cristoffel symbols 

_ fc _ 1 drndet( 7 ^) _ 1 ^ dcp a 
1 1 ~ 1 4fc " 2 9^ ~ 2 ^ 0?' ^ j 

a 

and their spatial derivatives 

^ = iV^ (2 33) 

a 

Tj do not depend on partial derivatives of angular variables 0j. Expressions for spatial deriva- 
tives of Cristoffel symbols are 



ax a ox K ox n ox K ox n 

(3) dc^dcf^ (4)p _^r_^,(5) r 0^.00 



(2.34) 



where 



a t a r> ij:rskn t „ t 1 ij; rs kn Q „ Q t ! 

^ ^-Tjj- r fc n — ~^(-A r i t A. r i5 n j -\- A r kA. r j5 n i 6^ ^ A an ^4 a ^,^4 r .jyl r jf S j , (2.35) 



^ r'yjrsfcn — r) -'AfcAsnAri-'^Tj'^ ^ q A a i c A cm A r iA r j5 sr e^' ^ , (2.36) 



ij;rkn — 2 ifc ^-£*r;ai T 6 ^-ac-^bi^rfic J "jn 

+^4* (s r;aj + e^-^ A^^-Br-fc) <5 in (2-37) 



2^afc e< ^ b ^"^an {^Br\bi-^-bj ~\~ -^-bi^r;bj J ! 
^ ^r'yjj.g^ = —(B r - a ] t B s - a i + A a f;.Esr;ai ~\~ Qai;rska]$jn~\~ 



2 

q ( B r - a fcB s\aj ~\~ A a fcE sr\aj ~\~ Qaj;rska j ^in ~p:(Qij;rskn + Qji \rsknj 

where 

Qij;rskn — ^ b ^ ^B J .. a j ; .A tl7 j.B S j; ) j.Aj ) j + A a kB T - an B g-foAbj-^ 
AakA an E S r\bi-Abj ~\~ A a fcA an B s-biBr;bj I > 



(2.38) 



(2.39) 



and 



^ ^^ij;rskn — ^ i^r.sk-^-si ~i~ A s kB r - s i^j 5j n -\- — ^B r ; S kA s j + A s kB r - s j^ b{ n 
— e"^ ^ i^B r - a }^A an + ^4aA;^r;an^ ^sjr 

+ ^~^A an A ab B r . sb (A si 5 kj + A si fe) (2.40) 



2 
1 

2 
1 

2 

Three-dimensional Ricci tensor 



-^ak-^-an^ ^ ( y B r - s iA s j + j4 s jB r - S j) 



^ = -d-dJ + F ^- r ^% (2-41) 

written in terms of new variables will contain second order quasi-linear terms proportional to 
dx k tx n an< ^ dx k dx n P^ US l° wer order non-linear terms made of products of first derivatives Inl- 
and 5%. We can write it as 

ox K 

r> _(1) R d 2( ^V | (2)p <9 2 6> r , 

i.7 r Hj;rkn n u n „ > -^ijirkn n u n „ ' 

ax K ax n ox K ox n 42) 

(3) p d(j) r d(j) s | (4)p 89^86^ {5) dcf^de^ 

OX K OX n OX K OX n OX n OX K 

Using formulas 1)2,32)1 - 1)2. 40|) . coefficients in front of second-order terms can be written as 

^ ^ ^ij;rkn = 2 ^-^-rfc-^ri^nj "I - ^rfc^rj'^rai < ^ >a A- an A a kA r iA T j^ "^fc^i™^"'' (2.43) 

and ^ 

^ ^ -^ij-rskn — ~^-A r iA r j ^A s ) ( .A sn e t ^ r ^ s 5 sr j4 a /j A ( j W e < ^ r < ^' a '^. (2.44) 



Low-order term coefficients ( 3 > 4 ' 5 )i? can be trivially found from 1)2.32)1 - ()2.4U)1 as well. 



2.4 Trace and trace-free part of the first differential 

ESR representation of the metric 7^ leads to a natural decomposition of its first differential 
djij onto trace and trace-free parts. Trace-free part depends only on d9i whereas trace part 
depends only on dfa. 

We denote trace of a matrix in a tangential space as 

tr(aij) = Sijdij = an, (2-45) 

and trace of a tensor in curved space as 

Tr{ aij ) = j ij aij . (2.46) 
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From (|2.2U|) we have 

tr(djij) = A ni C k;nm A mi d9 k + A ni e K d(f> n A ni = C k . nn d9 k + 5„„e*"#n = e K d<t> n . (2.47) 
Therefore, 

djij = d~i\j + d 7 */, d^j = A ni e^"d<j) n A nj , dr^ = A ni C k]nm A mj d6 k , (2.48) 
where superscripts t and tf denote trace- and trace- free parts of djij, respectively. Analogously, 

Tridjij) = Y j djij = Akie'^Akj (A ni C k . nm A mj dO k + A ni e^ n A nj dcj) n ] 

n n 

so that in a curved space we have as well 

d lij = d 1 f j + d 1 T j F , d7? = ^o-, d 7 f/ = d^. (2.50) 

3 Einstein's equations in ESR form 

Our starting point is a standard ADM 3+1 formulation which consists of an evolution part 



(2.49) 



dt 
dt 

and constraints 
where 

and we defined 



2aKij + Vifc + Vjfo 
ViVja + a (Rij + Idj) + (3 m V m Kij + K m V 3 f3 m + K jm V^ 

C = R + K 2 - K tj K ij = 0, 
C % = VjKf - ViK = 0, 



(3.1) 



(3.2) 



R = 1 ij R ij , K] = ^ m K jm , K = K™, (3.3) 
Kij = KKij - 2K?K jn . (3.4) 



3.1 Evolution part and constraints 

To keep up with the first time - second space derivative form of ()3.1|) . we introduce additional 
ESR variables, rates of deformation ijj k and angular velocities uj k , 

d(j>k , de k 

Vk, -ttt = w fc . (3.5) 



dt r ^ dt 

Our goal is to rewrite ()3.1|) in terms of 4> k , 6 k , tp k and Lo k . 

Taking time derivative of the first equation in (|3.1|) and combining it with the second 
equation, we get 



d 2 Hj _ _ 9 d" 



~ : - - - + 2aViVja - 2a 2 {//, ; + £«) + S, r (3.6) 
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where we defined 

S H = §- t ( V ^i + V i&) " 2a (/? m V m i^- + KinSfi™ + i^V./T) . (3.7) 

Next step is to rotate 1)3.6(1 into a coordinate system where 7^ is diagonal. Making use of ()2.21)l 
we obtain 

2a Aj n ylj m /C nm + 1~lij. nrn ip n ip m + T~Lij. nm Lo n uj m -\- l~L^j. nm ib n 9 m -\- j4j m ^4j n 5 mn . 

(3.8) 

First term on the left-hand side of 1)3.8)1 has non-zero elements only on the main diagonal 
whereas the second term has only off-diagonal non-zero elements. We will use this to obtain 
separate equations for tpi and Ui, but first we need to derive explicit expression for the right- 
hand side of 1)3.8)1 in terms of new variables. 

From the first equation in 1)3.1)1 . making use of 1)2.10)) we obtain 

Kij = -—Aaie^Aaj^a - —A ai C k . ab A bj u k + — (V^- + Vjft) , (3.9) 



so that 



1 . . , 1 . „ . 7 ifc 



iff = ^ k K ki = -—A kl A kj ip k - —e-^A ai C k , ab A b3 uj k , +±- (V^ + V fc ft) , (3-10) 
2a 2a 2a 

and 

if = ^ = -^y>fc + — Vi/3 fc = + ^-VA (3.11) 
2a a a a 

fc 

where 

^ = 1^V*- (3-12) 

k 

From ()3.9j) we see that first term on the right-hand side of 1)3.8)1 is simply 

— 2"t^" Ai n Aj m K nm = — ^6ije^ >t ipi + Ck-ijtoi^j , (3.13) 

and 1)3.8)1 thus can be rewritten as 

'ch/'i dm a \ / 9cjfc <91na 

2<(xAi n Aj m X7 n X7 m Q! 2a Ai n Aj m R nm -\- Ai m Aj n S mn 

2a Ai n Aj m JC nm -\- T~Lij. nm lj) n ll} m ~\~ 'Hij-nm^n^m ^■ij-.nm^Pn^m- 

Using 1)53)1 - (|5~TT)) . we obtain 

lb / 

( e ^ C ^'fc;cfe^6iCr;ce^.ej l - i; r'- t; A; + 6^" AdA c jll) c l/j c + ^ciCy-jce^j^rV'c + C ktCb A b iA c j1p c UJ k ^ . 

(3.15) 



(3.14) 
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Rotation of /C^ gives 

2a 2 A ri A sj fCij = 8 ra e^ r (ipijj r - ip 2 ) - C n . rs u n Nj r +i/j s -ifA - e~^ c C fc;cr C n;cs a;„u; fc . (3.16) 

First term in (|3.16|) has only diagonal, second only off-diagonal, and the third has both diagonal 
and off-diagonal non-zero elements. 

Consider first a diagonal part of (|3.14l) . Gathering diagonal contributions from the last four 
terms on the right-hand side of (|3,14|) we obtain a set of evolution equations for deformation 

rates ipi, 



dt dt 
where 



Gathering off-diagonal contributions in (|3.14|) we obtain 



1 1 in ■ 



dip d\w (y. ( 1 \ 

ipi = 2a 2 e'^'Ai n A im i - V n V m a - R nm j - ipipi + ^i-,kn^n^k + A im A in S, 

(3.17) 

^i;kn = ^ c ~ ^* J Ai m Ai r Bk\cmBn\cri ^ ] ^ ' i\kn = 0. (3.18) 



C k ;ij [ g^UJ k + lpUJ k 



2(xAi n Aj m \7 n \7 m Q, lot Ai n Aj m R nm -\- H,^j. nm ui n io m -\- ^l^j. nm ip n uj m -\- Ai m Aj n S mn , 

(3.19) 



where 



O 1 . —A., A. f> ,R ( p <t>i+<t>i-<t>c _ J>c 
iL ij;nm — ■ n it- n js ±J n;ct J - J m;cs V c c 



(3.20) 



Bji'isB m -j s (^e J +6 8 ) 6 l E mn -ikAjk a J E mn -jkAikj 

^ij;nm = ■^ik-^m;jk^ " ' ($in &jn) AjkB^ikG^ 1 (5jn &in\ (^-21) 

O 1 . = O 1 - O? = O 2 - O 2 = D C\ 22") 

* L ij;nm * b ji;nmi * L ij\nm * L jt;nmi "i!;nm 

Finally, solving equations (|3.19j) we obtain a set of equations for angular velocities uji, 
du>k din a 

~al ^~ u k + W^k = 

dt dt 

^kr \^rij \ i^(^Ai n Aj m V n V m O Iol Ai n Aj m R nm + ^l^j. nm Ul n LO m + Qij.^ nm 1p n U} m + Ai m Aj n S mn ^j . 

(3.23) 

Equations ()3.17j) . (|3.23|) together with two equations (|3.5|) constitute an evolution part of an 
ADM system transformed to ESR variables. 

Rewriting constraint equations is straightforward but the resulting formulas are rather 
complicated and we do not present them here. To derive these formulas, one has to substitute 
expressions for the Riemann tensor (|2.42|) , extrinsic curvature ()3.9j) , p.lUj) , (|3.11j) , and Cristoffel 
symbols (|2.31|) into (|3.2|) . In rewriting the momentum constraints, one has also to differentiate 
Kj and K with respect to x k , which in turn requires differentiation of terms containing rotation 
matrices and their first derivatives. Expressions of second-order derivatives of Aij are given 
in TTM . 
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For the energy constraint, one can also obtain a simple formula expressing Co through 
right-hand sides of the evolution equations for ipi. Summing up (|3,17|) and using ()3.11|) we 
obtain 

°* 91n % = -a 2 fi? + ^-Vv n V m ^ + 



dt dt V a / (3.24) 

2^ nm V n p m - ( 7 nm V n /3 m ) 2 + \p nm S nm . 

First and second terms in brackets in ()3.24j) can be replaced using the expression (|3.2j) for the 
energy constraint, 



— ~ M — r = ~ a \Co+K nm K - —7 V n V m a ] + 



dt dt 



2^7 nm V n /? m - ( 7 nm V„/3 m ) 2 + \^ nm S n 



(3.25) 



The latter equation can be inverted to obtain the energy constraint as 

■z r ly^^fc , 9\na 2 f 1 \ 

« = — - 2 ^ H Qt~y ~ \ K nm& ~a V n V m a I + 



2^ 7 nm V„/? m - ( 7 nm V n /? m ) 2 + 



(3.26) 



We see that in ESR formulation the energy constraint depends linearly on time derivatives of 
ipi. Setting Co = in ()3.26j) will define a plane in the three-dimensional space of in which 
the evolution of a constrained solution must occur, regardless of a choice of gauge. The unit 
normal to the plane is a vector with components 77^ 73^' ^ distance from the plane is 
given by the actual value of — a 2 Co- Modifying (|3.17|) by adding the term — a 2 Co/\/3 to each of 
the right-hand sides will project numerical evolution onto a local Co = plane, and may help 
with improving the accuracy of the evolution along a (generally curved) surface of constrained 
solutions. 

We also note that differentiation of ()3.11|) gives 

dK 1 fd-ijj , dlncA d fj nm V n p m \ . w , 

so that we can rewrite (j3.24[) . assuming Co = 0, in a more familiar form as 

^ = aK nm K nm - 7 nm V n V m a + "shift - dependent terms". (3.28) 
3.2 Case of degenerate eigenvalues 

As was mentioned in the beginning of the paper, decomposition (|2.1|) is always possible but it 
is not unique if 7ij has degenerate eigenvalues. In case of two degenerate eigenvalues, rotation 
angle in the direction orthogonal to the corresponding eigenvectors is arbitrary. In case of 
three degenerate eigenvalues all three rotation angles are arbitrary. For example, a flat three- 
dimensional hypersurface with Cartesian coordinate system on it is a case of triply-degenerate 
eigenvalues. In all of these cases, the matrix Cjj in (|2.19j) becomes degenerate and cannot be 
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inverted. Equation (|3.23|) then cannot be used to determine Below we describe how 
should be determined when eigenvalues are degenerate. 

Consider first the case of triple degeneracy, fa = fa = fa- Using (|2.14|) we can rewrite 
(t2~T3l as 

C k;ij = A in B k . jn (e^ - e*). (3.29) 

We see that in a triply-degenerate case all three Ck-ij = 0, and (|3.19|) becomes a system of 
three algebraic equations 

2o:Ai n Aj m \7 n \/ m oi lex Ai n Aj m R nm -\- 17^- .^ nm uj n uj m -\- ^^j. nm ijj n u m -\- Ai m Aj n S mn — 

(3.30) 

with respect to Ui. We assume that conditions for lapse a and shift are given, and that fa, 
9i, ipi, and U{ are known on a three-dimensional hypersurface as a result of previous evolution. 
Thus, first, second, and last terms in (j3,3()|) are known, and £lh. nm and £lh. nm which depend 
on fa and 9i are known as well. Then ()3.3(J|) is a set of three quadratic equations for Ui which 
are automatically satisfied at a degenerate point by continuity of u>i. Differentiation of (|3.30|) 
with respect to t will give us expressions for ^ which will depend on fa, 6i, ipi, Ui and on 
The latter values are given by (|3.17j) . Thus we can uniquely determine ^ as functions of fa, 
6i, ipi, u>i, and their spatial derivatives, which was our original purpose. 

The case of two degenerate eigenvalues can be treated similarly. Without loss of generality, 
assume that (pi = fa. Then all elements of symmetric matrices Ck-ij are zer ° excluding Ck-iz = 
Cfc;3i 7^ and C k -23 = Cfe;32 ^ 0. Equations (|3.19|) then reduce to one algebraic equation with 
respect to 



laAinAimS nS mCL ~ 2c? Ai n A2 m R n m + 



^Uinm^n^m + ^12;nmV ; n u; m + Ai n A2 m S nm — 0, 

plus two partial differential equations 



(3.31) 



-^r 77— + V'Wfc ) =2aA ln A 3m V n V m a - 2a 2 A ln A 3m R nm + 

dt dt ) (3.32) 

^l^nm^n^m + ^13;ramV'n'- L 'm + A\ n A 3m S nm , 

(dco c^ln a \ 
-kt k— + ipu k I = 2aA 2n A 3m V n V m oi - 2a 2 A 2n A 3m R nm + 
dt dt J (3.33) 

^23;nm^n^m + ^23;nmVVik-'m + A2 n A 3m S nm . 

Differentiation of (|3.31|) allows us to determine one time derivative function of 

other variables and two other derivatives, Substitution of into (|3.32|) and (|3.33|) 

will then give us two equations to explicitly determine -^f. Using them we can find ^§f-- 
A similar approach should be used in case of degenerate eigenvalues to determine SRE 
variables and their derivatives from 7^ and its derivatives on the initial hypersurface. We see 
from the above consideration that starting from degenerate initial conditions with arbitrary 
rotation angles 9i does not mean that rotation velocities tOi and accelerations can also be 
set arbitrary. Angular velocities must be determined from (|3.3Uj) . and they will depend among 
other things on lapse and shift conditions that we wish to use for the evolution. 
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4 Conclusions 



In this paper we formulated a tensorial exponential transformation of a three-dimensional 
metric of space-like hypersurfaces embedded in a four-dimensional space time in terms of expo- 
nential stretch-rotation, or ESR variables (Section 2). We derived formulas relating derivatives 
of the metric with the corresponding derivatives of these variables. A 3+1 system of Einstein's 
equations, formulated in terms of ESR variables, describes time evolution of the metric at ev- 
ery point of a hypersurface as a continuous stretch and rotation of a local cartesian coordinate 
system in a tangential space (Section 3). 

We want to mention that certain exponential transformation of variables were used in 
GR before for special symmetric cases. Such transformations sometime simplify analytical 
operations, and final results may sometime be formulated in a more compact form The 
transformation (|2.1|) . (|2.4|) . ()2.8j) considered in this paper is the most general exponential 
transformation which can be carried out for a three-dimensional metric without assuming any 
symmetries. 

A potential advantage of the ESR formulation is a change it introduces to the structure of 
right-hand sides of Einstein's equations by eliminating contravariant four-dimensional metric 
g ab as a multiplicative term in these equations. For a scalar non-linear hyperbolic equation 
considered in the removal of a term g -1 , where g is a scalar equivalent of the metric, lead 
to a dramatic improvement in long-term stability of numerical integration. Consider right-hand 
sides of (|3.17|) and (j3,23|) . Ricci tensor Rij and other terms written using ESR variables consist 
of parts that either do not have multipliers oc e~^ i or have multipliers of the type oc e^ -9 ^. 
That is, they either do not have j 13 multipliers or have multiplier that are ratios of different 
eigenvalues of 7^. In front of the right-hand sides There are multipliers of the type a 2 e _< ^ 
in front of the right-hand sides. We must recall that the lapse is a g 00 component of the 
four-dimensional metric g a b and thus having this multiplier is equivalent of having ratios of the 
zero-th (negative) eigenvalue to other three (positive) eigenvalues of a four-dimensional metric. 
We see that the net result of transformation is replacement of g ab with ratios of four-metric 
eigenvalues. We hope that the right-hand sides of the evolution equations (|3.17|) and (|3.23f) 
may behave better than the right-hand sides of the original ADM equations. 

We want to discuss now differences between the ESR formulation of Einstein's equations and 
other versions of 3+1 formulations of GR. Zelmanov [Zj, in his version of a 3+1 decomposition, 
chooses a congruence of time- like lines (motion of local reference frames), and then at each 
event on a time-line draws a small element of a three-dimensional space orthogonal to that line 
(a local three-dimensional reference frame). Einstein's equations in his formulation describe 
physical accelerations, deformations, and rotation of these reference frames (and physics of 
matter if space is not empty). There are no global three-dimensional slices embedded in a 
four-dimensional space-time in his formulation. In a classical ADM 3+1 formulation of 
GR and subsequent 3+1 formulations, a four-dimensional space-time is split into a family of 
three-dimensional space-like slices and a universal "time" is introduced to label these slices. A 
global coordinate system exists on each slice, and a three-dimensional metric of slices evolves 
with time according to Einstein's equations and pre-determined gauge conditions. In tetrad 
formulations of GR [HJ, a family of orthogonal 4- vectors is introduced and Einstein's equations 
are formulated in terms of variables projected onto tetrad components. 

The ESR formulation presented here is closest to a standard 3+1 ADM formulation. How- 
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ever, instead of metric and its derivatives, it uses variables which can be thought of in the 
following way. A three-dimensional metric defines an orthogonal triad of space-like vectors di- 
rected along the main axis of the metric tensor at each point of a hypersurface. ESR variables 
4>i and 9i tell how a three-dimensional coordinate system on a hypersurface must be rotated and 
re-scaled at each point so that it will locally coincide with a three-dimensional Cartesian co- 
ordinate system defined by the orthogonal triad at this point. The evolution of ESR variables 
is equivalent to the evolution of the three-dimensional metric and is governed by Einstein's 
equations written in terms of these variables plus the pre-defined gauge conditions. There is 
no orthogonal tetrad in this formulation. 

From the above discussion it is clear that ESR 3+1 formulation presented here can be 
modified and extended by introducing new variables such as spatial derivatives of the metric, 
and by addition of various combination of constraints, similar to modifications introduced to 
a standard ADM formulation before [2]. Hyperbolicity and stability properties of ESR system 
and its modifications will be studied in subsequent publications. We plan to use this new 
system for numerical solutions of certain problems of GR. 
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